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ABSTRACT 

We use narrowband imaging (FWHM - 70 A) to select a sample of emission line galaxies 
between 0.20 < z < 1.22 in two fields covering 0.5 sq. deg. We use spectroscopic follow-up 
to select a sub-sample of Ha emitting galaxies at z ~ 0.24 and determine the Ha luminosity 
function and star formation density at z ~ 0.24 for both of our fields. Corrections are made 
for imaging and spectroscopic incompleteness, extinction and interloper contamination on the 
basis of the spectroscopic data. When compared to each other, we find the field samples differ 
by Aa = 0.2 in faint end slope and Alog[L*(erg s"')] - 0.2 in luminosity. In the context 
of other recent surveys, our sample has comparable faint end slope, but a fainter L* turn- 
over. We conclude that systematic uncertainties and differences in selection criteria remain 
the dominant sources of uncertainty between Ha luminosity functions at this redshift. 

We also investigate average star formation rates as a function of local environment and 
find typical values consistent with the field densities that we probe, in agreement with previous 
results. However, we find tentative evidence for an increase in star formation rate with respect 
to the local density of star forming galaxies, consistent with the scenario that galaxy-galaxy 
interactions are triggers for bursts of star formation. 

Key words: surveys - galaxies; luminosity function, mass function - galaxies; starburst 



1 INTRODUCTION 

It is now widely accepted that the amount of star formation in Uni- 
verse as a whole has increased since the formation of the first galax- 
ies, peaking around re dshifts z ~ 2 - 3 and subsequently declining 
by a factor of ten (e.g. lHopkinsll2004l and references therein). Cos- 
mic star formation history provid es strong constraints on models of 
galax y formation and evolution ( |Peietal.llT999l :[ Somerville et al.l 
l200ll) . because it di rectly traces the accumulation of stellar mass 
and metal fraction jPei & Falllll995l ; iMadau et alj[T9 961 to their 
present-day values jCole et al.ll200ll : IPanter et al.ll2003b . Its rapid 



decl ine over the past 8 Gyr is consistent with "downsizing" scenar- 
ios ( ICowieetal.[|l996l) in which the more massive galaxies have 
produced their stell ar mass at earlier times than the less massive 
galax i es ( H eavens et al.l |2004| ; lluneau et al.l |2005| ; | Thomas et al.l 
I2OO5I ; iFardal et al.l 12006'). The star formation history of the uni- 
verse has also been used to constrain a l lowable stellar initial mass 
functions jBaldrv & Glazebrookll2003l : iHopkins & Beacomll20o3) 



* Based on observations made with ESO Telescopes at the La Silla Obser- 
vatory (Programmes 67.A-0063, 68.A-0363 and 69.A-0314) and the Anglo- 
Australian Telescope, 
t Email: westra@mso.anu.edu.au 
:j: Email: heath@aao.gov.au 



and c osmic supernova rates dOal- Yam & Maozll2004l : lDaigne et al.l 
I2OO6I) . 

Star forming galaxies exhibit a strong UV continuum courtesy 
of newly formed OB stars in sites of star formation. This newborn 
population can be infer red from the UV directly (e.g. iTrever et al.l 
ll998l ; lLillvetai]|l996l) or through a host o f indirect calibrators 
spread across the electromagnetic spe ctrum jRosa-Gonzalez et al.l 
l2002l ; ICondoi]|l992l ; ISchaereillioOol) . At low redshifts the most 
direct calibrator - and of the optical calibrators the least affected 
by internal extinction - is the Via recombinatio n line, whi ch emits 
when stimulated by ionising UV radiation (e.g. iKennicutl|[T998.) . 

Narrowband surveys at optical wavelengths have long been 
recognised as a powerful way of yielding large samples of emis- 
sion line_^alaxies, including those selected by Ha at redshifts z ^ 
0.4 jLvet al. 2007; Pascual et al. 2007; lones & Bland-Hawthor^ 
1 200 lb . Thev are advantageous in that they select galaxies in exactly 
the same quantity that they seek to measure, and are optimised for 
the detection of the faint emission line signatures indicative of star 
formation. Narrowband surveys also have the advantage of a sim- 
plified selection function, with filters that probe only a very nar- 
row redshift slice, thereby yielding a volume limited sample at a 



Lya at high redshift ( 


Aiiki et al.l2003l:lHu et alj2004l;lRhoads et al.l 


|2004 iGawiser et al.l 


20061). as well as Ha, H/3. TOml and [Oiil at 
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lower redshifts jpuiita et al.ll2003l : iHippelein et alj|2003l ; ILv et alj 
l2007h . 

Here we describe a survey for Ka emission line galaxies at 
z ~ 0.24, found as a by-product of the Wid e Field Lyman Alpha 
Search (WFILAS: Iwestra et alJlioO^ . l2006h . The resulting sam- 
ple has been utilised to determine the Ha luminosity function at 
z ~ 0.24 and its associated co-moving star formation density. In 
Section|2]we describe the selection of candidates using narrow- and 
broadband imaging. In Section[3]we detail follow-up spectroscopy 
used to identify the nature of the emission and test completeness of 
the sample. In Section|4]we derive the Ho- luminosity function for 
galaxies at z ~ 0.24 and explore its variation with the local environ- 
ment in Section |5] A summary and concluding remarks are made 
in Section|6] 

Throughout this paper we assume a flat Universe with 
(n^.^A) = (0.3,0.7) and a Hubble constant Hq 
70kms~' Mpc~'. A ll quoted magnitudes are in the AB system 
dOke&Gumill 1983111 . 



2 CANDIDATE SELECTION 
2.1 Narrowband imaging 

The observations were done with the Wide Field Imager (WFI) on 
the ESO/MPI 2.2 m telescope at the Cerro La Silla Observatory, 
Chile. The WFI consists of a four by two array of 2k x 4k CCDs 
giving a total field size of 34" x 33"with pixel scale of 0'.'238 per 
pixel. Imaging dat a were taken from the Wi de Field Lyman Alpha 
Search (WFILAS : rWestra et alj|200^.l2006h , a wide-field narrow- 
band survey des igned to find Lyman-tt emitte rs at z ~ 5.7. We re- 
fer the reader to I Westra et al.l (j2006, hereafter IPaper J) for a more 
detailed description, but give the important features of the survey 
below. 

Three fields spaced around the sky were observed in three nar- 
rowband filters (FWHM = 7 nm) centred at 8 10, 8 17 and 824 nm, an 
intermediate width filter {FWHM = 22 ma) centred at 815 nm and 
broadbands B and R. For one of the fields with missing 817 nm data 
it was not possible to apply the selection criteria uniformly and so it 
was excluded from this analysis. The two fields used were the well- 
studi ed Chandra Deep Field South (CDFS; e.g. 'Rosati et al."2002"; 
iRix et al. 2004) and the COMBO- 17 Sll field (Wolf etal. 2003). 
The width of our narrowband filters is essentia lly half that of other 
surveys (e.g. iFuiita et ai] |2003l :l Lv et al. I l2007h with a correspond- 
ing reduction in background and enhancement in the contrast of 
observations of emission line galaxies. Table [T] gives an overview 
of the emission lines redshifted into these narrowband filters, the 
associated luminosity distances and co-moving volumes. 

The data were processed using a combination of standard 
IRAljl routines (mscred) and some custom designed for our data. 
Image frames were bias-subtracted, flat-fielded and background- 
subtracted. A fringe pattern present in the intermediate band and 
narrowband images, which remained after the flat-fielding, was re- 
moved using a fringe frame created from 10-30 science frames. Fi- 
nally, an astrometric correction was applied using the USNO CCD 

' "'AS = -2.5 log /,, - 48.590, where m^g is the AB magnitude and /,, is the 
flux density in ergs s"' cm"^ Hz"' 

IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 



Astrograph Catalogue 2 (UCAC2: Izacharias et"al]|2004l) and the 
IRAF-task msccmatch with a resulting RMS of < 0".' 15. 

To ensure the quality of the final deep images we only in- 
cluded frames with a seeing of less than 5 pixels (= V!2) and without 
significant fringing. The images were weighted according to their 
exposure time and combined using the IRAF mscstack routine re- 
jecting deviant pixels. 



2.2 Photometry and completeness corrections 

We used SExtractor (version 2.3.2; iBertin & Amoutsll996l) in dou- 
ble image mode to do create the initial source catalogues. Each 
resulting catalogue contains the photometry for the sources in all 
6 filters. Sources were selected when at least 5 pixels were 0.8 cr 
above the noise level in the narrowband image used for detection. 
All photometry was measured in apertures with a 10 pixel diameter 
(= 2'.' 4). Paper I describes the procedure in detail. 

Detection completeness was determined using galaxy number- 
counts in each of the narrowband images as a function of AB- 
magni tude and that of the H ubble Deep Field (HDF) in the FilAW 
filter j Williams et alj[l996h . Completeness is defined in this in- 
stance as the ratio of the number of detected galaxies to that of 
expected, and the completeness correction is its reciprocal. The ex- 
pected number counts were fit by a simple linear function over the 
magnitude range [20, 25]. For all the objects that are selected as our 
candidates this correction is less than 0.1 %. 



2.3 Selection criteria and star/galaxy disambiguation 

The following four criteria were applied to select our candidate 
emission-line galaxies from the initial source catalogues: 

(i) the narrowband image used as the detection image must have 
the most flux of all the narrowband images and the source must 
have a 4 cr detection or better in the detected narrowband; 

(ii) there must be at least a 2cr detection in the intermediate 
band image; 

(iii) the broadband image R needs to have a 2(T detection or 
better; 

(iv) the emission line flux calculated from the narrowband im- 
ages should be Fii„e 5^ 10"'^ erg s"' cm"^. 

Criterion |(iii)| removes emission line objects with a very low 
continuum. These objects were cl assified a s Lyo- emitters at z ~ 5.7. 
The Lyo- emitters are discussed in lPaper j . The emission line fluxes 
that we use in this paper were measured from the narrowband 
photometry. The background (or underlying continuum) was de- 
termined by averaging the flux measured in the two narrowband 
images that were not used for the detection of the source. This was 
subtracted from the flux measured in the narrowband detection im- 
age, which is emission line and continuum flux combined. An aper- 
ture correction was calculated according to: 

C = max(0.2,erf|^jxerf|i^J) (1) 

and applied to the line fluxes. Here, C is the fraction of light of the 
object contained within the 10 pixel aperture, a and b are the profile 
width along the major and minor axes, respectively (assuming that 
the galaxy profile is adequately represented by a two-dimensional 
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Emission line 








Hot 


w 


[Oiii] 


[On] 


[Sii] 


Redshift range in A'gio 


0.229-0.239 


0.659 - 0.673 


0.610-0.624 


1.163-1.182 


0.199-0.210 


Redshift range in A'8i7 


0.239-0.250 


0.673 - 0.687 


0.624-0.638 


1.182-1.201 


0.210-0.220 


Redshift range in N^i^ 


0.250-0.261 


0.687 - 0.702 


0.638-0.652 


1.201 - 1.219 


0.220-0.230 


Dl (Mpc) 


1203.1 


4081.9 


3726.3 


8158.4 


1045.1 


VCDFS (10^ Mpc3) 


9.4 


60.6 


53.6 


137.2 


7.3 


Vsil (lO^Mpc^) 


8.3 


53.2 


47.1 


120.6 


6.4 



Table 1. Redshift coverage, luminosity distance D^, and co-moving volume for each emission line in each of our narrowband filters A'gio, Nm and ^824 using 
the central wavelength and FWHM of each filter (70 A). The CDFS and S 1 1 fields span differing volumes ( VcDFS and Vs 1 1 > respectively). For [Oiii] we used 
the wavelength of the [Oiii] A 5007 line and for [On] and [S n] the average wavelength of the individual lines of each doublet. 



Gaussian) and erf (x) is the error functioiQ. To ensure that the fluxes 
of certain large objects were not over-corrected, we limited C to 
at least 0.2. Dividing the calculated emission flux by C gives the 
emission line flux Fn^^ used in criterion [(iv)[ 

The emission line flux limit of 10"'* erg s"' cm"^ is a fac- 
tor of two higher than the detection limit of our earlier search 
for high redshift Lya emitting ga laxies u sing the same imaging 
data (Fiinii, = 5 x 10"" erg s"' cm~^: |PaDer t) . This is because we are 
no longer limited by the night-sky background, but rather by the 
brightness of the object continua. This limit was chosen in part 
to ensure that emission line candidates were within the sensitiv- 
ity limits of our follow-up confirmation spectroscopy. We note that 
we used a flux limit rather than an equivalent width cut-oif. The 
lowest equivalent width values as determined from the naiTowband 
imaging in our candidate sample is ~ 1 A, with a peak at ~ 3 A. 

Stars represent a significant fraction of contaminants. We 
found that standard star/galaxy classification from SExtractor 
works satisfactorily for objects brighter than R = 21. However, it 
breaks down for the large number of faint (R > 21) objects. There- 
fore, additional criteria were applied. We examined the size of the 
objects (major and minor axes), in combination with their shape 
(the ratio of the major and minor axes) as additional star/galaxy 
discriminants. Since this size/shape information could potentially 
lead to the unwanted removal of unresolved line emitting galax- 
ies, we used an additional cut in (B - R) colour as a safeguard to 
prevent thifl We decided to restrict the size/shape discrimination 
to sources with (B - R)> 1.4 based on the (B - R) colour distri- 
bution of Ha emitters at z ~ 0.24 and [Sii] emitters at z ~ 0.21 
obtained from spectroscopic observations. This size/shape/colour 
criterion was added after initial spectroscopic follow-up to im- 
prove removal of stellar contaminants. Furthermore, the colour of 
(B - R) = 1.4 corresponds to the model of an instantaneous star- 
burst with a n age of ~ 1 Gyr. We d etermined this colour using 
GALAXEV teruzual & Charlotll2003h . We finalised our stellar se- 
lection criteria as follows: 

(i) the SExtractor CLASS.STAR parameter is > 0.95 and < 21. At 
R>21, sources are too faint for SExtractor to reliably distinguish 
between stars and galaxies; 

(ii) the SExtractor A_iMAGE and b_image parameters (the profile in 
pixels along the major and minor axes, respectively) are < 4 pixels, 



^ The error function is defined as erf (x) = fe '^dt 



* We note that this colour selection will also reject QSOs. This is of no 
consequence to our selection of a star forming sample. 
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Figure 1. Distribution of observed (B-R) colour for naiTowband candidates 
satisfying various criteria for stellarity: (a) SExtractor class.star > 0.95 and 
i? < 21 (forward cross-hatching), (b) SExtractor < 1.06 and (B - 

R) Ss 1.4 (backward cross-hatching) and (t) bright stars showing dift'raction 
spikes or ghost reflections (horizontal cross-hatching). The thick histogram 
shows the combined distribution of Ha and [S 11] galaxies (z ~ 0.24 and 
z ~ 0.21, respectively) from our full emission line sample subsequently 
through follow-up spectroscopy (Section [3.U , scaled by 0.2. 



the ratio of these parameters is g-™^™ < 1.06 and the object has 
a(B - R) colour > 1.4. This is redder than almost all star forming 
galaxies at z ~ 0.24; 

(iii) the object showed obvious imaging artefacts, such as 
diffraction spikes or ghost images, in any of its thumbnails. 

Figure [T] shows the distribution of spectroscopically observed ob- 
jects that satisfy these criteria as a function of observed (B - 
R) colour for the CDFS field. The forward cross-hatched his- 
tograms represent objects satisfying criterion [(1)] the backward 
cross-hatched those for criterion |(ii)[ and the horizontal cross- 
hatched those for criterion [(iii)[ The histogram outlined by the thick 
solid line represents the observed (B - R) colour distribution of se- 
curely confirmed Ha and [Sn] emitters (z ~ 0.24 and z ~ 0.21, 
respectively), by way of comparison. All objects selected in this 
way were deemed to be stellar and removed from the candidate 
list. Finally, all candidates were inspected to remove sources that 
were contaminated by image artefacts. 

From initial candidate numbers of 786 and 848 for the CDFS 
and Sll fields respectively, 414 and 513 candidates were removed 
because they met one or more of the stellar criteria. Our final sam- 
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Figure 2. Observed (B - R) colour distribution for various sets of emission- 
line galaxies within our sample: (a) Ha at z ~ 0.24 {forward cross- 
hatching), (b) [Sii] at z ~ 0.21 (backward cross-hatching) and (c) single- 
line emitters of indeterminate origin (horizontal cross-hatching). 



pie yielded 372 candidate emission-line galaxies for the CDFS field 
and 335 for the SI 1 field. 



3 SPECTROSCOPIC FOLLOW-UP 
3.1 Observations and reduction 

The emission-line selection criteria established in Section 12.31 
are sensitive to almost any galaxy with emission lines that have 
been redshifted into the wavelength range of our narrowband fil- 
ters, and are bright enough to be detected. The one exception is 
Lya, which does not yield detectable flux blueward of the Ly- 
man limit and hence in our broadband images. The main emission 
lines to expect in our narrowband filters are (from bluest to red- 
dest), [On] ^^3726,3728, Hy ^4342 (although usually too faint, or 
too much underlying absorption), H/3 /14863, [Oiii] /1/14959,5007, 
Hot /16564 and [Sii] /i/i6733,6718. Since the goal of this paper is 
to establish the star formation density at z ~ 0.24, we concentrated 
only on those gal axies detected as Ha. Alternative approaches by 
other groups (e.g. lLvetal]|2007h have separated objects based on 
their broadband colours. Unfortunately, in the case of [S ii] galax- 
ies (z ~ 0.21) the colours are indistinguishable from those with Ho' 
iz ~ 0.24) due to their similar redshifts. Figure|2]shows how the Ha 
and [S n] galaxies occupy the same range of colour [(B - R) > 0.5] 
given their near-identical redshifts. Based on this, we classify all of 
the single-line emitters outside this range [(B - i?) < 0.5] as likely 
[On] line-emitters at z ~ 1.2. It is worth pointing out that when the 
[S n] doublet falls inside our narrowband filter set, an extra volume 
of about 50 % of the volume probed by Ha can be explored. Unfor- 
tunately, the fluxes of [S ii] and Ha are not sufficiently correlated 
to p ermit star formation density determinations from the [S ii] line 
(e.g. lKewlev et"Ii]|200lh . and so it was not used. 

Our approach was to target as large a sample as possible of our 
candidates to test how successful our candidate selection was. An 
additional aim was to measure the fraction of the observed candi- 
dates with Ha in our narrowband filters. To do this, we ensured 
that the spectroscopic sample was representative of the narrow- 
band sample as a whole. A two-sided Kolmogorov-Smirnov test 
yielded probability levels of 99.8 % and 49.3 % for the CDFS and 



S 1 1 fields, respectively. Once measured, we applied the determined 
fraction to our entire sample of candidates in each field. 

The spectro scopic data were taken with AAOmega 
dSharp et al.l |2006|) . an optical multi-object spectrograph. It is 
fibre-fed from the prime focu s of the Anglo-Au stralian Telescope 
(AAT) by the 2dF facility jLewis et alll2002h to a dual-beam 
spectrograph, which in our case was used with spectral ranges 
3800-5700 A and 5700-8700 A. The resolving power was 
SA = 3.5 A in the blue arm and 6A = 5.3 A in the red arm. It has 
392 fibres available to observe spectra of objects within a 2 degree 
field of view. The fibres have a minimum placement separation 
of 30", although the actual limiting separation depends on the 
orientation of fibre buttons when placed on the field plates. For 
fields with a high density of targets, such as our 0.5° x 0.5° fields, 
only ~ 250 fibres could be allocated per configuration, due to such 
placement limitations. In general, the number of fibres allocated 
depends upon the target distribution in the field and the choice of 
algorithm in the configur^f] so ftware. We found tha t using the 
Simulated Annealing algorithm jMiszalski et al. 1 1200^ allowed a 
larger fraction of fibres to be allocated to candidates than the older 
Oxford algorithm. 

The data were taken during four separate runs. The first obser- 
vations were done in classical mode during 2 nights, 2006 March 
23 and 24. During this run the Sll field was observed. The other 
three occasions were done in service mode on 2006 October 10, 
2006 November 10, and 2007 March 26. During these runs both 
fields were targeted. We used the 580V and the 385R volume phase 
holographic (VPH) gratings for the blue and red arm, respectively. 
Table [2] summarises the observations. In total, 301 and 255 candi- 
dates were observed in the CDFS and Sll fields, respectively. 

Basic spectral reductions, including bias-subtraction, flat- 
fielding and wavelength calibration were done using the 2dF re- 
duction pipeline drcontrol^. The final one-dimensional spectrum 
for each object was obtained by averaging the reduced spectra of 
the object in the different observations using our own IDL scripts. 

The spectra of several stan dard stars (LTT 7379, LTT 7987 
and CD-32 9927; iBessellll 19991) were taken during the final night 
of the 2006 March run and were reduced in the same fashion as the 
science data. System throughput as a function of wavelength was 
derived using each standard star and its sensitivity curve. These 
curves were scaled to a common level and averaged to give the 
overall sensitivity. This was applied to all the science spectra to 
flux calibrate each relative to one another. Unfortunately, absolute 
flux calibrations are very difficult to do reliably with fibre-based 
spectrographs, due to the changing configurations of the fibres and 
the effect this has on their throughput. For this reason, we used the 
line ffuxes measured from our narrowband photometry rather than 
the fibre spectroscopy. 

3.2 Spectroscopic completeness 

We used a Monte-Carlo simulation that combined the background 
of real spectra of our securely confirmed Ha emitting galaxies with 
transplanted and scaled emission lines to asses our spectroscopic 
completeness as a function of line ffux. We took the spectrum of 
each Ha emitter and fitted the Ha and [Nii] lines together with 
the continuum. Each line was fitted by a Gaussian and the galaxy 

^ configure and drcontrol are software packages produced and 
maintained by the AAO. These packages can be obtained from 
|ftp : //ftp ■ aao . gov . au/pub/2df| 
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observing 


field 


number of 


tot. exp. 




dates 


observed 


configurations 


time (sec) 


seeing (") 


2006/03/23 


Sll 


2 


11,700 


0.9-1.5 


2006/03/24 


Sll 


3 


14,400 


1.3-1.8 


2006/10/10 


CDFS 


1 


9,900 


1.8-2.2 


2006/11/10 


CDFS 


1 


11,700 


1.2-1.5 


2007/03/26 


Sll 


1 


6,300 


2.5 



Table 2. Details of the spectroscopic follow-up observations. 



continuum (or background sky) was approximated by a first order 
polynomial. The line centres were parameterised by redshift. The 
widths of the [N ii] lines were set equal and t he flux ratio bet ween 
the red and blue [N n] lines was fixed to 2.96 liMendoziJll983h . The 
remaining fit parameters were left unconstrained. The model of the 
Ha- [Nil] complex was subtracted from our data, leaving only the 
underlying noise. To the noise, we added a randomly scaled ver- 
sion of our model with a random offset in wavelength. We then at- 
tempted to re-identify any emission line. We did this multiple times 
for each secure Ha emitting galaxy. 

This exercise demonstrated that it was possible to identify at 
least 90 % of the galaxies at a line flux of log Fiine = - 16.0 (Fii„e in 
ergs"' cm"^) for all spectroscopic runs. In Figure[3]the recovered 
fraction as a function of line flux is shown for the CDFS and Sll 
fields. The uncertainties indicated in Figure [3] were derived using 
the following relation: 



o"f, 



VA'to,(A?d=, + 2)(A'to,-/Vto + 1) 



(2) 



MotCMo, + 3) 

where crf^ac is the calculated uncertainty, /Vtot the total number of 
objects in that bin and N^a is the number of o bjects which have a 
detection of the emission line (after Eq. 4 from ljones et alj|2006h . 
The spectroscopic completion rate as indicated in Figure [3] is well 
fit by a function of the form 



77(F) 



exp[-y(F - F^f°] F < F, 
1 F^ Fc 



(3) 



where y represents the speed at which the function drops off and F^ 
is the flux at which the function reaches 1.0. 



3.3 Ha emission line fraction 

In almost all cases the spectra of confirmed emission-line galax- 
ies should show additional emission lines elsewhere except cases 
of Lya at z ~ 5.7 (which are filtered out through their absence of 
B and R flux) or [On] at z ~ 1.2. This is demonstrated by Fig- 
ure |4l where we show the stacked spectrum of all our confirmed 
Ha and [S n] galaxies in the CDFS. Ha is usually accompanied by 
the [Nn] /i^6550,6585 and [Sii] /i^6733,6718 doublets, whereas 
Hj8 and the [Om] /1/14959,5007 doublet are almost always seen to- 
gether. Our spectral resolution {R ~ 1500 at 8150 A) is not enough 
to fully resolve the [On] doublet, but high enough to show it as 
broader than a single emission line. The signal-to-noise ratio of the 
spectra is not always high enough to clearly determine if a line is 
broad (in this sense) or not. Alternatively, these galaxies could be 
Ha emitting galaxies with all other emission lines too faint to be 
detected. 

There are a few galaxies which show only one emission line. 
Although we expect many of them to be [On] emitters at z ~ 1 .2, we 
cannot rule out the possibility of single-line Ha galaxies at z ~ 0.24 




-16,5 -16.0 
bg Flux (ergs s^i cm-^) ^ 



Figure 3. Spectroscopic completeness as a function of line flux for the 
CDFS (top) and Sll (bottom) fields as derived from a Monte-Carlo sim- 
ulation of artificially generated emission lines. See text for details. 



in without the use of additional information. In Figure |2] we show 
the observed (6 - R) colour distribution of galaxies in the CDFS 
where the emission line in the narrowband filters has been con- 
firmed as Ha or [S ii] through the presence of additional lines. We 
also indicate the colour distribution of galaxies for which we have 
only one emission line feature. Some of the single-line detections 
are bluer than the combined Ha/[S ii] distribution. We therefore 
identify all single-line galaxies with (S - /?) < 0.5 to be [On] emit- 
ters at z ~ 1.2 and those with (B - R)> 0.5 to be Ha emitters at 
z ~ 0.24. 

Of the candidates for which we have spectroscopically con- 
firmed an emission line (189 and 117 out of the total 301 and 255 
observed in the CDFS and Sll fields, respectively), just under a 
half are Ha at z ~ 0.24, a quarter are [Sn] at z ~ 0.21, roughly a 
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Figure 4. Mean spectrum of emission-line galaxies from the CDFS field. Spectra from 114 galaxies between z = 0.i9 and 0.27 were de-redshifted before 
stacking. The most prominent features have been labelled. The apparent absorption feature just bluewards of the [Oi] /!6302 line is the remnant of the telluric 
A-band of the individual spectra being de-redshifted and stacked. This spectrum was used to fit the emission Hnes to derive the mean extinction as described 
in Section [T4l Only red arm data from AAOmega (observed wavelength ~ 5700 - 8700 A) are shown. 



sixth are H/? or [Oin] at z ~ 0.6 - 0.7 and the remainder are [On] at 
z~ 1.2. 

Figure [5] shows the fraction of confirmed Ha emitters in our 
full spectroscopic sample as a function of narrowband flux. It peaks 

around logFyne 15.3 (with F],„^ in ergs"' cm"^), below which 

increasing numbers of [On] galaxies at z ~ 1.2 begin to dominate 
the counts. Each point in Figure [5] has a minimum of 10 galaxies 
per bin and a minimum binwidth of 0. 1 dex. The uncertainties in 
the Ha fraction per bin have been calculated using Eq. |2l where 
A'det now represents the number of galaxies with confirmed Ha. We 
fit a Gaussian of the form 



f (F) = ax exp 



20-2 



+ h. 



(4) 



where F,. is the flux central to the peak, cr and a are its width and 
height, and is a zero-point offset. The resulting fits are shown in 
Figure[5] 

We decided to fit both the CDFS and SI 1 fields individually, 
given the likely differences between the field samples due to cos- 
mic variance. Given the relatively narrow range of volume probed 
through each emission line, we expect over- and underdensities 
at the different redshift intervals to change the relative numbers 
of galaxies as a fu nction of flux jjones & Bland-Hawthonj|200ll : 
IPascual et all200lh . 



3.4 Extinction corrections 



Star forming regions are some of the dustiest galaxy environments, 
making correction for internal obscuration necessary. Many emis- 
sion line s urveys apply a general extinction correction of Ann- ~ 1 
(e.g. iTress e & Maddox 1998; Fujitaetal. 2003). However, it has 
been shown t hat there are large variations in ex tinction between 
alaxi es (e.g. Ijansen et al] 1200 ll) . Furthermore, iMassarotti et al.l 
200 ih state that applying an average extinction correction always 
underestimates the true extinction correction. Since our spectra 
cover a large wavelength range (3800-5700 A in the blue and 
5700-8700 A in the red) we are able to observe Ha and H/3 simulta- 
neously. We therefore calculate the extinction individually for each 
galaxy through Ha and Hy6 when both lines are detectable. The 
signal-to-noise ratio is not always high enough to show H^S clearly 




i5,5 -15.0 
Loq flux (ergs s"' cm"^) ^ 



-16.0 -'5,5 -15.0 -14.5 

Log flux (ergs a-' cni"^) ^ 

Figure 5. The Ho- fraction of our candidates for the CDFS (top) and S 11 
(bottom) fields. The data have been binned to have a minimum of 10 galax- 
ies per bin and a minimum width of 0. 1 dex. The dotted line is the four- 
parameter fit to the data points. The horizontal error-bars represent the width 
of the bins and the vertical bars the uncertainty in the fraction calculated us- 
ing Eq. Q. 
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Figure 6. Colour excess E(B-V) and Ahu (assuming k{iia) = 2.47) as a 
function of B-magnitude as determined by measuring the Balmer decre- 
ment of the averaged spectra in the CDFS. The crosses indicate the colour 
excess without using any correction for stellar absorption in Ho' and HyS. 
The diamon ds indicate the c olour excess using the same stellar absorption 
correction as lHopkins et"ai] {2003) with EW(Ha) = 1.3 and EW(H8) = 1.6. 
The dotted and dashed line are the linear fits to the respective points. 



in emission. Therefore, we grouped available spectra according to 
the S-magnitude of the source, obtained an average spectrum, and 
measured the Balmer decrement value from these. 

The colour excess E(B - V) can be calculated using 



E{B - V) ■ 



2.5 log Ral3 

KH/j) - k(Ha) 



(5) 



where Ra/s is the ratio of the observed value of the Balmer decre- 
ment to its theoretical value, and A:(H/3) - k(lla) is the differen- 
tial extinction between the wavelengths of H/3 and Ho-. The the- 
oretical value for the Balmer decrem ent is 2.87 (fo r T = 10'* K 
and case B recomb ination; Table 2 of lCalzettHbOOll . which uses a 
ICardelli et al.l 19891 extinction law) and the value for the differential 
extinction is 1.163. This assumes k(V) = 3.1 and k(lla) = 2.468. We 
adopt these values throughout the rest of this paper. 

In Figure |6] we plot the resulting values for E(B - V) as a 
function of the S-magnitude for two cases: without and with cor- 
rection for absorption due to the underlying stellar population. The 
AAOmega spectra have a resolution of ~ 5.3 A throughout the red 
arm meaning that we are unable to resolve the HjS absorption line 
directly. If we assume no stellar absorption, the colour excess has 
values up to E(B - V) ~ 1, corresponding to Ann- = 2.5 mag (or 
Av = 3.1mag) using A^ = k(A) x E(B - V). This is far higher 
than the average extinction of Aua ~ 1 as assumed elsewhere (e.g. 
iTresse & Maddoxlll998l : lFuiita et al.ll2003h . If we instead adopt the 
median equivalent widths for stellar absorption in Ha and H/3 as 
measured by Hopkins et al. ( 2003), 1.3 and 1 .6 A respectively, then 
the average extinction as shown in Figure [6] is roughly Ahb ~ 0.85. 
We note that there is a trend of a decreasing extinction with increas- 
ing apparent magnitude (see Figure[6l(. Observe that our sample has 
a restricted range in redshift, making apparent magnitude B a proxy 
for absolute magnitu de Mg. Sim ilar trends of change in E(B - V) 
have been found by I Jansen et al.l j200ll) . We attribute this trend to 
the fact that either fainter (and therefore smaller) galaxies poten- 
tially contain less dust, or the H/3 flux might be overestimated in the 
mean spectrum of the faintest galaxies as a result of a low signal-to- 



noise ratio of the H^ line. We derive an extinction of Ahu = 0.9(0 
from the Balmer ratio in the mean spectrum of all emission-line 
galaxies as shown in Figure |4] Since the trend might be due to a 
low signal-to-noise ratio of the H/J line, we use a constant value 
throughout to correct for extinction. 



4 LUMINOSITY FUNCTION AND STAR FORMATION 
DENSITY 

4.1 Derivation and fit 

With the final emission line catalogue in hand, and the various se- 
lection and completeness effects accounted for, our approach to 
calculating the Ho- luminosity function is as follows. We take our 
measured distribution of line emitters (all emission lines from all 
redshifts) from the narrowband candidate sample and apply the 
spectroscopically measured fraction of Ho- emitters as a function of 
flux (Section[33}. We correct for incompleteness in both the spec- 
troscopic identifications (Section [3.2| ( as well as the original nar- 
rowband imaging. The corrections for the latter are less than 0.1 % 
(Section [2. 2t . Finally, we correct our line fluxes for the effects of 
extinction (Section r3.4| (. 

Figure [7] shows separate luminosity functi ons for both the 
CDFS and S 1 1 fields. We fit a Schechter function l lSchechteill973) 
to the data points using a minimised fit. The Schechter function 
is given by 



0(L).L = ^-(ilf exp(-il).(il) 



(6) 



where 0' represents the normalisation constant of the galaxy den- 
sity, a the faint end slope, and L* the characteristic luminosity 
where the Schechter function rapidly declines at bright luminosi- 
ties. We used a Levenberg-Marquardt method for finding the min- 
imum x~ fit to the binned data-points, courtesy of the IDL routine 
mpfitfun from the Markwardjj library. Since the three parameters 
a, L* and ^* are highly correlated, we used the correlation matrix 
and the partial derivatives of the Schechter function to calculate the 
formal uncertainty in the integrated luminosity density X, 



d£d£_ 

dxj dxj 



(V) 



Here, X| , X2 and corresp ond to the Schechter parameters a, log L* 
and log(^* l lCowaniri998l) . Vjj is the covariance matrix, which re- 
lates to the correlation matrix p,j as V,y = PijcrjCTj. CTj is the formal 
uncertainty in the parameter. We list the resulting values of the 
parameters and the formal uncertainties, together with the correla- 
tion matrices in Table[3] 

The luminosity density over luminosities L ^ Liin, can be cal- 
culated by integrating Eq.|6l yielding 

£ = 0*LT(Q' + 2, j^) . (8) 

In the case where limiting luminosity Ly,^ = 0, the luminosity den- 
sity reduces to X = (^*LT(q' + 2). Using the Schechter parame- 
ters and uncertainties given in Table [3] with logLyn, = 40.6 (Lum 



* Alternatively, the extinction law of ICalzetti et al.l )2OO0l) gives Ahu = 
1.18, where k{V) = 4.05, Jt(Ha) = 3.325 and kiHfi) - k{Ua) = 1.163. This 
correction results in Ha fluxes roughly 20 % higher than the values used in 
the text. 

^ Maintained by C. Markwardt at |http://cow.physics.wisc.edu/~craigin7] 
idl/idl.html. 
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Figure 7. Top: luminosity function for Ha galaxies at z ~ 0.24 for the CDFS (left) and Sll (right) fields. The solid Hne in each of these panels is the fit to 
the data points, while the dotted line indicates the fit of other field for reference. Bottom: Confidence levels for the parameters a, L' and 0* of the CDFS (left) 
and Sll (right) fields. Contours are drawn for each plane in which one of the parameters is held constant. The 1, 2 and 3 cr contours indicated correspond to 
68.3%, 95.4% and 99.7% confidence limits, respectively. 



CDFS 

a log L* log <p' 

-1.33 ±0.34 41.43 ±0.22 -2.23 ± 0.32 
' 1.00000 -0.91020 0.96458 ~ 
-0.91020 1.00000 -0.97268 
0.96458 -0.97268 1.00000 



sn 

a log L* log <p' 

-1.11 ±0.51 41.24 ±0.25 -2.28 ± 0.33 

' 1.00000 -0.90948 0.95099 " 

-0.90948 1.00000 -0.96826 

0.95099 -0.96826 1.00000 



Table 3. Schechter parameters for the Ha luminosity functions for each field determined using a Levenberg-Marquardt minimisation. The correlation 
matrices p,j for each are shown below. 



in ergs ', corresponding to our survey flux limit) gives logX = 
39.17+°;?^ and 38.86!jJ;|i in ergs"' for the CDFS and Sll fields, 
respectively. The uncertainties are calculated using the correla- 
tion matrices in Table [3] If we instead use the Ha luminosities of 
the galaxies directly and sum over all, we obtain 39.22;^j5[J^ and 
38.86;^g[j3 for CDFS and Sll, respectively. The uncertainties in 
this case are the square-root of the sum in quadrature of individual 
galaxy luminosity uncertainties and does not take into account Ho- 
emission line fraction uncertainties and, as such, are lower limits. 



4.2 Comparison to previous surveys 

In Figure[8]we compare our Schechter fits to the results of other sur- 
veys using Ha as a measure for star formation. The survey parame- 
ters are summarised in Table|4] We restricted the comparison to Ha 
surveys with z g 0.40 in order to limit the systematic uncertainties 
which play into the comparison when different star formation in- 
dicators are involved. It can be seen that there is a large range in 
each of the Schechter parameters between surveys, a ranges from 

1.1 to -1.6, logL* from ~ 41.3 to 42.2 and log0* from -3.7 

to -2.2. Some of these surveys cover different redshifts to those 
in our survey. The wide span of the parameters could be attributed 
by evolu t ion of the luminosity func tion, as has been suggested by 
iHopkinj j2004h and lLv et al.l ( l2007h . who compare surveys over a 
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Figure 8. Comparison of the Schechter functions derived for our two fields 
(thick solid lines) and those of other surveys. The Schechter function of each 
survey has been drawn over the luminos i ty rang e where data was available. 
The other surveys are Tresse & Maddox 1 1998, dash). Sullivan et al.| j200d 
dash dot),[Fuiita et al. (2003, long dash), Hippelein et al. (2003, thick shoil 
dash), and lLv et al.l <2007l z = 0.24 sohd; z = 0.40 thick long dash). The 
individual Schechter parameters are given in Table|4] 

wider redshift range using different indicators. However, a number 
of systematic uncertainties exist between surveys that could also 
attribute to the scatter between the luminosity functions. We now 
explore each in turn. 

The details of galaxy selection inevitably vary from survey 
to survey. For example, Tresse & Maddox ( 1998) have selected 
their galaxies from an /-band selected sample, while iSullivan et al.l 
( I2OOQI used UV imaging to select theirs. It is well known 
that galaxy selection based on different passbands results in a 
different faint end slope of the galaxy luminosity distribution 
(Madgwick et al. 2002; Jones et al. 2004). Passbands that favour 
bluer and/or star forming galaxies generally yield higher faint end 
counts and thus steeper slopes. This undoubtedly has a similar in- 
fluence on the faint end slope of the Ha luminosity function. 

It is also important to note that any survey using an equivalent 
wid th selection (or eq uiv alently, a narro wband-broadband colour, 
e.g. lFuiita et"ai]|2003l and lLv et alj|2007l) . unlike our survey which 
rather applies a flux limit, will tend to be biased against galaxies 
with low equivalent widths. This will affect mostly the selection of 
galaxies with a high star formation rate per unit continuum, such as 
early type spirals ( Kennicutt 1992) and galaxies with low Ho- flux in 
general. The Ha luminosity function, of course, only characterises 
the line flux on its own. 

In Section [J!4l we discussed the amount of extinction correc- 
tion for our survey and concluded that it agrees with values found 
by other surveys. However, there is still a large spread in the ex- 
tinction values. A range of = 0.5 - 1 .8 is typical of those found 
dLv et alj I2OO7I : iKennicuttI 1 199 8!. and references therein), which 
translates directly into an uncertainty of 0.3 in logL*. The excep- 
tion is when all galaxies have individually been corrected for ex- 
tinction, which imposes large observational overheads. None of the 
surveys indicated in Figure[8]have been able to do so. 

Some surveys have only a limited spectroscopic follow up on 
their candidates, or none at all (Fujita et al. 2003; Ly et al. 2007). 
Both of these surveys use additional colour criteria to distinguish 
between Ha an d other line emitting galaxies at other redshifts. 
ILv et al.l ( I2OO7I) estimate that there is about 50 % contamination 
of [Oiii] galaxies into the Ha sample o f iFuiita et alj j2003h based 



on empirical colour selection using spectra from the Ha waii Hub- 
ble De ep Field-North. Spectroscopy on several sources in lLv et al] 
( l2007h shows that slight contamination of higher redshift emission 
line galaxies occurs in their Ha sample. In our own sample, as 
we noted in Section [3711 there is a large sample of galaxies that 
has been selected on their [S 11] lines, which would have otherwise 
been mistaken for low redshift Ha had we relied on colour selec- 
tion on its own. Furthermore, the fraction of contamination by other 
emission line galaxies varies significantly with obs erved line flux 
jjones & Bland-Hawthornl2001 ; |Pascual et alj200ll) . Hence, spec- 
troscopic observations of all or a large representative sample of the 
candidates is vital in understanding the amount of contamination 
by galaxies at different redshifts. 

Spectroscopic observations also allow flux corrections for the 
[Nil] /1/16550,6585 lines, which straddle Ha with an observed sep- 
aration of ~ 44 A at z ~ 0.24. A more detailed analysis of [Nn] 
is described in Section [4.4. 11 Fluxes quoted in this paper do not 
include this correction, unless otherwise stated. 

Cosmic variance has widely been cited as a major contri butor 
to the differences between various surveys (e.g. lLv et al.l2007h . We 
are well-placed to test the impact of this given that we have ob- 
served two distinct fields that have been subjected to identical se- 
lection and analysis. We have estimated the contribution of cosmic 
variance to the mean object densities given by the luminosity func- 
tions i n Figure [8] Following the prescription of ISomerville et al] 
( l2004h we determined the relative cosmic variance al for several 
Ha surveys. The estimate of cr,, is an upper limit as our survey 
has the shape of an elongated prism, wh ile the derivation is for a 
spherical volume jSomerville et al.|[2004l) . The cosmic variance is 
calculated by = baxm, where b is the bias parameter (defined as 
the ratio of the root variance of the halos and the dark matter) and 
(Tq^ the va riance of the d ark matter. Using a number density of 
0.05 Mpc"^ dLv et al.ll2007l) yields a bias ofb ~ 0.7 for all surveys 
at z < 0.40. The corresponding variance over our survey volumes 
(of 9.4 X 10^ and 8.3 x 10' Mpc^) is a^u ~ 0.7 and thus cr,, = 0.49. 
This translates to an uncertainty in log0(L) of -I-0.2/-0.3, which is 
ample to account for the difference between the luminosity func- 
tions of the two fields. 

Many of the narrowband surveys exhibit similar uncertainties 
which are sufficiently large to account for the differences between 
each other. Table |5] shows resulting uncertainty in the number den- 
sity due to the cosmic variance for a sample of narrowband sur- 
veys with we l l-defin ed survey volumes. Despite the low redshift, 
lGallegoetalJ ( ll995h span a large enough volume that their uncer- 
tainty due to cosmic variance is somewhat lower than the surveys 
at higher redshift. Comparing the uncertainties Alog^(L) to the 
spread of luminosity functions in Figure[8l we observe that cosmic 
variance is one of the dominating factors in the determination of an 
average Ha luminosity function at these redshifts. 

Finally, we make the observation that there is a high degree of 
correlation between the three Schechter parameters. This is clearly 
demonstrated by the confidence limit contours in the bottom panels 
of Figure [7] and the correlation matrices in Table|3] 



4.3 Star formation density 

The amount of extinction-corrected Ha luminosity from an Hii re- 
gion is directly proportional to the quantity of UV ionising flux 
produced by newborn stars. As such, it can be used to estimate the 
number of new stars and hence the star formation rate. We can thus 
derive global star formation densities from the Ha luminosity den- 
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Reference/field 



Redshift 



logL* 



log <l>' 



Oallego et al. (1995) 


0.022 


± 0.022 


-1.30 ± 


20 


41 


87 


± 


08 


-2 


79 


± 


20 


Tresse & Maddox (1998) 


0.200 


±0.100 


-1.35 ± 


06 


41 


92 


± 


13 


-2 


56 


± 


09 


Sullivan et al. (2000) 


0.150 


±0.150 


-1.62 ±0 


10 


42 


42 


± 


14 


-3 


55 


±0 


20 


Fuiita et al. (2003) 


0.242 


± 0.009 


-1.53 ± 


15 


41 


95 


± 


25 


-2 


62 


±0 


34 


Hiopelein et al. (2003) 


0.245 


± 0.022 


-1.35 






41 


.45 






-2 


.32 




Perez-Gonzalez et al. (2003) 


0.025 


± 0.025 


-1.20 ±0 


20 


42 


43 


± 


17 


-3 


00 


± 


20 


Lv et al. (2007) 


0.0735 ± 0.0075 


, 0.0855 ± 0.0055 


-1.59 ±0 


02 


42 


05 


± 


07 


-3 


14 


± 


09 


Lyet al. (2007) 


0.242 


± 0.009 


-1.71 ± 


08 


42 


20 


± 1 


24 


-3 


70 


± 1 


06 


L^etal. (2007) 


0.401 


±0.010 


-1.34 ±0 


06 


41 


93 


± 


19 


-2 


75 


± 


16 


This paper, CDFS 


0.245 


±0.016 


-1.33 ± 


34 


41 


43 


± 


22 


-2 


23 


± 


32 


This paper, S 11 


0.245 


± 0.016 


-1.11 ± 


51 


41 


24 


± 


25 


-2 


28 


± 


33 



Table 4. Values for the parameters of the Schechter functions shown in Figure[8] V is in erg s ' and ip* in Mpc ^ . After Table 5 in lLvetaljj2007h . 



Reference 


Redshift range 


Sky area 


Co-moving volume 


tr,. 


Alog<^(L) 






(sq. deg.) 


(10^ Mpc-3) 




(0(L) in Mpc-3) 


pallego etal. (1995) 


z ^ 0.045 


471 


3.3 X 10^ 


0.21 


-I-0.1/-0.1 


Vujita et al. (2003) 


0.233 < z sS 0.251 


0.255 


3.9 


0.56 


+0.2/-0.4 


HiDpelein et al. (2003) 


0.238 < z < 0.252 


0.086 


1.4 


0.70 


+0.2/-0.5 


Lv et al. (2007) 


0.233 < z < 0.251 


0.255 


4.7 


0.63 


+0.2/-0.4 


This paper, CDFS field 


0.229 < z < 0.261 


0.262 


9.4 


0.49 


+0.2/-0.3 


This paper, S 1 1 field 


0.229 < z < 0.261 


0.230 


8.3 


0.49 


-I-0.2/-0.3 



Table 5. The survey geometries for a sample of narrowband surveys with well-defined survey vo lumes alongside their ro ot cosmic variance and the associated 
uncertainty in the number density. Root cosmic v ariance was calc ulated using the prescription of lSomerville et all j2004h assuming bias h = 0.7 and a number 
density of line emitters of 0.05 Mpc"^^, following lLv et al.l )2007l) . 



sities of Section l4ni 
lKennicuttl ( fT993) , 



We use the star formation rate calibration of 



p (Mo yr"') = 7.9 x lO^^^LCHa) (ergs s"' cm"^) , 



(9) 



which assumes a Salpeter initial-mass function, case B recombina- 
tion and an electron temperature of 10* K. 

In the case of some surveys jGallego et all I2OO2I : 

iHi ppelein et al ] I2OO3I : IPerez-Gonzalez et alj booi " this paper), 
the faint-end slope of the Ha luminosity function is poorly 
constrained, thus having important consequences for the integrated 
luminosity density. To overcome these, and in order to make a 
fair comparison, we calculate the star formation density of other 
Ho- emission line surveys at the same redshift by assuming a 
common fixed limit rather than integrating from zero luminosity. 
We choose = 0.33 Mq yr"', which corresponds to the limit 
of our survey (logFum = -16.0 with Fu^ in ergs^'cm"^, or 
logLiiin = 40.6 with Lii,„ in ergs"'), and avoids faint-end extrapo- 
lations or assumed faint-end fits of some other surveys. This yields 
logp(L > Liim) = -2.24!;;;]^ and -1.93!;j;°g for the S 11 and CDFS 
fields, respectively. 

Our two fields are indicated in Figure |9] The other results in- 
cluded in this Figure are derived in the same way as described with 
Eq. l[8j in Section [43] We included only star formation densities 
from surveys based on emission lines and transformed onto the 
same cosmology. The majority of the se points were calculated us- 
ing the compilation of ILv et aU ( |2007|). We also incl uded the least- 
squares fit to the z < 1 points of Hopkins! ( |2004|) as a point of 
reference. Note that this fit assumes Liin, = 0. 



Observe that the star formation density in both our fields 
agrees quite well with other Ha emission line surveys at the same 
redshift. Nevertheless, there is a difference of almost 1 dex between 
the highest and lowest v alue for the star fo rmation density. The 
highest value comes from lFuiita et alj l2003h . which (according to 
ILv et al.ll2007l) suffers from contamination of higher redshift emis- 
sion line galaxies, pushing their value upwards accordingly. Ob- 
serve in Fi gure [9] that we h ave also plotted the star formation den- 
sity fits of lHopkinsI ( l2004h which, unlike the points, make use of 
star formation density values integrated down to zero luminosity. 
This serves to illustrate the extent to which extrapolation of the 
faint end fit affects the final determination of star formation density: 

typically up to < 50 % for a 1.3 (larger for steeper values). As 

discussed earlier, the luminosity functions of several surveys have 
ill-constrained faint end values. 

Obviously the same systematic uncertainties discussed in Sec- 
tion |4j2] will also play a role here. Furthermore, since we compare 
the star formation density over a larger redshift range, other emis- 
sion line star formation indicators have been used (most notably 
[On]), thereby introducing their own sources of systematic uncer- 
tainty In the case of [On], extinction corrections are larger and 
its star for mation rate calibrato r depends on the abundance of the 
ionised gas l lKewlev et alj2004h . Corrections for both can be made 
with spectra covering Ha, Hj8 and [Om], as well as [On]. However, 
at redshifts z <; 0.5 these lines are progressively lost from the opti- 
cal, giving rise to uncertain ties of up to .4 in log(SFR), whe n ap- 
plying the lKennicutJ i ll 9921) calibrations jKewlev et alj|2004l) . Be- 

i ond this, emission line analyses are pushed into the near-infrared 
Glazebrook et al.ll999l : lDohertv et al.l2006h . where brighter night- 
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Figure 9. Star formation density as a function of look-back time derived 
from emission line surveys, where the Schechter function has been inte- 
grated from the star formation rate corresponding to the flux limit of our 
survey, 1 X 10"'* erg s"' cm"^ (0.33 Mq yr"' ). The solid symbols represent 
the star formation density derived from the Ha line, the open symbols from 
either the [On] or [Oiii] line. The solid circles are the star formation den- 
sity from the CDFS and S 1 1 fields of this pap er (top and bottom symbol, 
respectively). Other data are Fuiita et al.l ^20031 open and soli d diamonds), 

1 Sullivan et al.l i200(i open and soUd upward-poi nting triangle), [Tresse et al] 

solid downward-pointing trian gle), the ILv et al.l feOOTl open and 
sohd sq uares), [H ippelein et jfl j2003l open and solid righ t-pointing tri- 
angles ), iGallego et alj 1 199a solid left-pointing triangle), iGallego et alj 
i2002L open left-p ointing triangle), T resse & Maddox ( 199^ solid upward- 
pointing star) and IPerez-Gonzalez et al.l 12003 . solid downward-pointing 
star). The dotted and dashed line are the least-squares fit from Figures 1 and 

2 of / Hopkins 2004), respectively. They are not corrected for the fact that 
iHopkin s ( 2004) integrated the Schechter function down to L = erg s"' and 
are indicated for comparison purposes only. The parameters used to make 
this Figure aie given in Table |4] 

sky background and instrument thermal contributions increase the 
difficulty of making observations. 



4.4 Minor contaminating effects 

We now turn our attention to some additional sources of contami- 
nation for which we have not included any correction due to their 
minor nature. These are (i) the effect of the [N n] /1/16550,6585 lines 
on the measurement of Ho- flux, (ii) the presence of active galactic 
nuclei (AGNs) in the sample, and (iii) foreground [Sii] emitters. 
We consider each of these effects in turn. 



4.4.1 Na rrowband Ha flux and [Nil] 

The filters used in this survey are sufficiently narrow (FWHM= 
70 A) that the forbidden lines of [Nii] ^^6550,6585 can affect the 
measurement of Ho- flux through narrowband imaging. The prox- 
imity of lines complicates the estimate of Ho- line flux in two ways. 
In the case where Ho- is central to the filter, both [Nii] lines con- 
tribute to the narrowband-derived flux attributed to the Ho- flux. In 
the case where Ha is near the filter edge, then an [N ii] line is ei- 
ther lost to an adjacent filter (where it contributes to the measured 
continuum) or lost to the filters altogether Due to the finite width 
of emission lines and the steepness of the filters, this is a gradual 
transition. 
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Figure 10. Measured narrowband Ha flux in each WFILAS filter as a func- 
tion of redshift, expressed as a fraction of the true Ha flux. The thick sohd 
line shows the measured Ha flux in the absence of [N ii] and approximately 
traces the filter transmission curve. The shaded envelope shows the effect 
of the [N ii] line over the range Hff/[N ii]tot = 2.30 (dotted) to 4.66 (dashed), 
centred on 3.00 (thin solid line). The vertical dotted lines indicate where the 
filter transmission is 50 %. 



Ha 


Reference 


, [Nil] 


H(T( narrowband) 


LNiiJ,„, 




H(r(tiue) 
Nil! Ns24 


2.30 


Pascual et al. (2007) 


-0.49 


1.04 


1.04 1.14 


3.00 


this paper 


-0.60 


0.99 


0.99 1.08 


4.66 


Lv et al. (2007) 


-0.79 


0.93 


0.93 1.00 



Table 6. Ratio of Ha to combined [N ii] line flux for a range of values in the 
literature (as indicated). The third column shows the corresponding values 
of log of the ratio of the [N ii] /i6585 to Ha line flux. The last three columns 
indicate the mean ratio of the Ha flux measured in the narrowband filters to 
the true flux where the filter transmission is ^ 50 %. 



We have calculated the Ha flux as measured from the nar- 
rowband imaging and compared it to its true value over a range 
of galaxy redshifts, taldng into account our filter setup. This was 
done by taking the spectrum of one of our galaxies and fitting the 
Ho- line, as well as the [Nii]- and [Sii] doublets. We then con- 
volved the fit with each of our filter profilej^ and calculated the 
line flux from the filters in the same way as for the survey. A range 
of different values of HQ-/[ Nn]iot have been used, since this ratio 
depends o n meta llicity (e.g. IOsterbrock|[l98^ : iKewlev et alj[200ll : 
Kauffmann et al. 2003). We hav e used tw o extreme values (2.30 by 
f ascual et al. 2007 and 4.66 bv iLv et al.|[2007.) to reflect the wide 
range of metallicities found in these galaxies. A third value of 3.00 
was measured from a high quality spectrum used for emission line 
fitting. 

The results are shown in Figure[TO] Average values of the ratio 
of measured to true Ha flux are indicated in Table [6] The narrow- 
band Ha flux overestimates the true flux by about 10 % when av- 
eraged over all redshifts pertaining to a specific filter However, the 
ratio can peak around 40 % in the innermost 15 % of filter coverage. 
This peak corresponds to the specific case of an idealised square fil - 
ter containing all three lines as calculated bv lPascual et alj f2007h . 

^ We have used Butterworth curves for the transmission profiles. The curve 
is given by T{A) = |l + {iTp^Y } ■ where Ac is the central wavelength, h 
is the FWHM of the filter and n controls the steepness of the filter edges. 
We assumed to be 8100, 8170 and 8240 A for Wgio, Ngn and ^524, re- 
spectively, and FWHM = 70 A and n = 3. 
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This is a worse case scenario that only occurs rarely in practise. In 
the vast majority of cases the effect of [Nii] is moderated by the 
sloping edges of a real filter profile, or complete absence of [N ii] 
from the narrowband filters altogether. Given the 70 A width of our 
filters and the ~ 50A width of the Hq'/[Nii] group, the chance of 
having all three lines in the same filter is uncommon. 

As the overall effect of [Nii] is approximately 10% (corre- 
sponding to 0.04 in log L), we do not make any correction for it. 

4.4.2 AGN contribution 

The presence of an active nucleus in a galaxy can contribute Ho- line 
flux in addition to th at due to normal star formation. For example, 
iPascual et al. Il l200lh have found approximately 15 % of their lumi- 
nosity density to be due to galaxies identified as AGNs. We com- 
puted the fraction of Ha contribution due to AGN in our sample 
using the li ne diagnostic rela t ions a s determined by Kewley et al. 
( 1200 ll) and iKaufi^mann eTal] ( |2003|) . We selected galaxies from 
both fields where the fluxes of the emission lines Hj6, [Oiii] /15007 
and [Nil] /^6585 have been measured with a signal-to-noise ratio 
> 2. The line ratios and the line diagnostic re lations are indi cated 
in the Baldwin-Phillips-Terlevich (BPT; Baldwin etalj|l98lh dia- 
gram of Figure nn 

T wo of the galaxies lie above the extreme starburst demarca- 
tion of iKewlev et al.l ( 1200 ll, soUd line) and have log([Nii]/HQ') > 
-0.6, which classifies them as AGNs. A third galaxy, also with 
log([Nii]/HQ') > -0.6, lies below this demarcation , but ab ove the 
pure star formation boundary of iKauffhiann et al] ( l2003l dashed 
line), making this galaxy a composite case. 

If we consider the total Ha contribution from the two AGNs in 
it amounts to 5 % of the total Ha flux from this subsample. Overall, 
this AGN contribution would result in a decrease in the star forma- 
tion density of logp = 0.02. If we include the composite galaxy as 
a third AGN, the decrease is logp = 0.04. 

4.4.3 [S n] contribution 

Surveys that only apply colour criteria to determine the nature of 
the emission line are unable to distinguish between Ha emitters 
at z ~ 0.24 and [Sii] emitters at z ~ 0.21 (Section [3TT] and Fig- 
ure [2}. The previous calculations of the Ha fraction assume that 
it is possible to distinguish between Ha and [S ii] emitters. Since 
this can only be done with spectroscopy, we have also derived the 
Ha luminosity function for the case where [S ii] emitters are taken 
to be Ha emitters. This gives an indication of the impact of having 
sole reliance on colour criteria and no spectroscopic follow up. Fig- 
ure [T2] shows the difference between Ha luminosity functions for 
the CDFS and Sll fields where the [Sii] galaxies were assumed 
to be Ha. The Schechter parameters that belong to the alternative 
Schechter functions are given in Table [7] 

The star formation density as determined from these hypothet- 
ical [Sn]-as-Ha luminosity functions down to our survey limit are 



logp(L > Llim) : 



-2.22:^;^' and- 



i.sr 



for the SI 1 and CDFS 



fields, respectively. In the Sll field there were only a handful of 
galaxies (4) identified at z ~ 0.21, while the CDFS contained a 
significant number (30). Hence the results of the CDFS field are 
more significantly affected and an increase in the star formation 
density log p in this field of 0. 1 2 (about 30 %) can be seen. This ex- 
ercise demonstrates how a foreground overdensity of [S 11] emitters, 
if present, can significantly influence the Ha star formation density. 
For this reason, it is imperative to have at least some spectroscopic 
follow up to a narrowband survey to make such situations obvious. 
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Stor forming 
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Log {[Nll]/Hc<) ^ 

Figure 11. BPT diagram for 35 galaxies where the fluxe.s of the emis- 
sion lines H/S, [Oiii] A5Q01 and [Nii] /16585 have been measured with a 
signal-to-noise ratio ^ 2. Indi cated by the solid line is the extreme star- 
burst demarcation of lKewle~e t al. ( 2001) and the dash ed line the demar- 
cation of pure star formation of iKautfmann et al J )2003h . The vertical line 
is drawn at log([Nn]/HQ') = -0.6. The median uncertainty is indicated 
by the error bai's. The galaxies indicated by the open circles lie above the 
iKewlev et al . (2001) line and are most likely influenced by AGN activity. 
A third galaxy, marked by an open square, lies within the two demarca- 
tions and is a composite source. T he remaining gal axies, shown by the 
closed circles, lie below the Kaut fmann et al ] i2003ll relation and left of 
the log([Nn]/HQ') = -0.6 line and are pure star forming galaxies. 
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Figure 12. Differences between our derived Ha luminosity distribution 
(Figure [tJ and that assuming all [Sn] galaxies at z ~ 0.21 to be Ha. The 
values for the Schechter parameters to the luminosity diagram including the 
[S n] galaxies are indicated in Tablel?] 



5 ENVIRONMENTAL PROPERTIES 

The suppression of star formation rates at the centres of clus- 
ters has been well established both through direct observation 
(Lewis etal. 2002; Balogh et al. 1997, 1998; Kodama et al. 200ll 
|2004 iGomez et al. 2003), as wefl as a changing mix of mor- 
phological types (^Dressler.il98Q) . Such high density environments 
provide a range of dynamical mechanisms whereby galaxy en- 
counters rapidly strip gas from any potential star forming galax- 
ies (e.g. Couch et al. 2001, and references therein). Recent obser- 
vations have suggested a continuation of this trend across struc- 
tures at larger scales and lower density enhancements than clus- 
ters dGomez et al.ll2003l ; iGrav et al.ll2004h . Accordingly we exam- 
ine our two fields for evidence of star formation rates that are driven 
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-0.97085 
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Table 7. Schechter parameters for the Ha luminosity functions where the [S ii] galaxies at z ~ 0.21 are assumed to be Ha galaxies at z ~ 0.24. They were 
determined using a Levenberg-Marquardt minimisation. The correlation matrices p,j for each are shown below. 



by either the general galaxy environment, or alternatively, the local 
distribution of star forming galaxies. 

Usually, the amount of galaxy clustering is expressed as a 
function of projected density 



(10) 



where r„ (in Mpc) is the distance to the «th (usually n = 10) near- 
est neighbouring galaxy with Mg < -19. In cluster environments 
the star formation rate ha s been obse rved to be quenched at galaxy 
densities above 1 Mpc"^ jLewis et"a i. 2002; Gomez et al. 2003). 

In Figure [T3Ta) we show the fraction of galaxies with a star 
formation rate exceeding IMoyr"', as well as median and mean 
star formation rate per galaxy as a function of the projected den- 
sity of the general galaxy population. This uses data for all of the 
spectroscopically confirmed star forming galaxies at z ~ 0.24 for 
both of our fields combined. The indicated errorbars in the two top 
panels were determined using the jackknife estimatofl while in 
the bottom panel they are the standard deviation. We also show the 
25th and 75th percentile values for each bin in the middle panel. We 
determined the projected density by using the usual rig measure of 
the tenth-nearest star forming galaxy to each ordinary galaxy. Ordi- 
nary galaxies were taken f rom the photome tric redshift catalogues 
of the COMBO- 17 survey jWolf et al.l2"003l . K. Meisenheimer, priv. 
comm.) as galaxies with Bab < 22 (corresponding to Mg = -19) 
between 0.21 < z < 0.29. As the thickness of the redshift slice 
influences the value of projected density, we scale it using the dif- 
ference in the thickness of the redshift slice of our survey and the 
average thickness of the Scr cluster volumes (where <r i s the veloc- 
ity dispersion of the cluster) used in lLewis et a D|2002). 

Since we did not target any known clusters with our fields, we 
expect that there will be little or no evidence for star formation sup- 
pression in our fields. Typically, the projected density for galaxies 
within the virial radius of a cluster is ~ 4Mpc"^ and at the cen- 
tre of some rich clusters can be as high as lOMpc"^ ^Lewis et al.l 
I2OO2I) . Indeed, as Figure [T3T a) shows, there is negligible change 
in the star formation rate per unit density for the galaxies in both 
our fields (noting that the highest density point is affected by poor 
number statistics). Furthermore, we confirm levels of star forma- 
tion that are typical for the range of typical field gala xy densities 
probed by our data as fo und by previous surveys (e.g. lLewis et al.l 
l2002l : lG6mez et al .l2003h . Generally, the distribution of star forma- 
tion rates in a given density bin is rather asymmetric, making the 
median a more reliable measure than the mean. 

In Figure [T3lb) we show the same measures as for (a), but as 



The jackknife estimator is calculated as follows. Let p(,) = 
p(xi, . . . , Ar;+i, . . . , .v„) be the value of the statistic with one ele- 
ment X, removed, and define p(.) = (l/'O Pw- Then (tJ = (n - 
1)/" Z"-i (P(i") ~ P(-))^ is th e square of the jackknife estimate of standard 
error jEfron & Gonell983h . 
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Figure 13. (a) The mean and median star formation rate per galaxy (in 
Mq yr"' ) and the fraction of galaxies with a star formation rate > 1 Mq yr"' 
as a function of the projected density Z|o of ordinary galaxies (taken from 
the COMBO-17 survey). The errorbars in the two top panels are the jack- 
knife estimates of the standard error, while in the bottom panel they are 
standard deviations. The small crosses indicate the values for each individ- 
ual star forming galaxy. The thick lines in the middle panel indicate the 25th 
and 75th percentile for each bin. The numbers in the bottom panel indicate 
the number of galaxies included in each point. Some of the individual galax- 
ies have values outside the range of star formation rates plotted and hence 
are not indicated, (h) Same as {a), but as a function of star forming galaxy 
density ZJJ . 



a function of projected density of the spectroscopically confirmed 
star forming galaxies at 0.23 < z < 0.26. There are roughly one- 
third as many star forming galaxies as not, and so we redefine the 
projected density in terms of distance to the third-nearest galaxy, 
E3. As a consequence, and Sio span a similar range of density 
values. We observe in Figure [T3l b) that star formation per galaxy 
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Figure 14. Spatial distribution of the spectrally confirmed Ha galaxies in 
both our fields (solid circles). The size of the circles indicates the star forma- 
tion rate of the galaxy derived from the narrowband flux and the grey-scale 
the redshift. The black dots are galaxies that have not been spectroscopically 
confirmed yet and have a colour 0.5 < (B - i?) < 1.3, which corresponds to 
the colour interval of our confinned Ho: galaxies (Figure |2). 



increases with increasing density. Noting again that the highest den- 
sity bin is affected by poor number statistics. Although not con- 
clusive, this is consistent with galaxy evolution scenarios that see 
galax y-galaxy interactions as triggers for bursts of star formation 
Kionso et al. 2004; Perez et al. 2006). 

To examine the apparent relationship between star formation 
rate and projected density of star forming galaxies, we plot the spa- 
tial distribution of our spectrally confirmed Ha galaxies in Fig- 
ure [14] The size of the points indicates their star formation rate 
and their shade of grey the redshift. Probable (but unconfirmed) 
Ha candidates are also shown. These were selected on the basis 
of colour (0.5 < (B - R) < 1.3; see Figure [2} and having either 
indeterminate or non-existent spectra. 

The distribution of star forming galaxies in the CDFS field 
(Figure [T3b) suggests a tendency for grouping of the star form- 
ing galaxies. However, the eye is remarkably good at making out 
patterns in noisy distributions and thus we sh ould be cautious in 
these interpretations (e.g. p35 of Peebles 19931) . On the other hand, 
the distribution of star forming galaxies at z ~ 0.24 in the S 11 field 



(Figure[74t is apparently less structured than the CDFS. Because of 
this, we infer that the trend of increasing star formation with rising 
density of star forming galaxies is largely attributable to the data 
from the CDFS field. This contrast between the fields can also be 
seen in differences in the Ha space-densities given by the two lumi- 
nosity functions in Figure [7] As a consequence, the star formation 
density of the SI 1 field is lower than that of the CDFS (Figure (Sj. 
This is due to the lower Ha fraction in the S 1 1 field compared to 
the CDFS field (Figure|5}, to the extend that can be seen given the 
more limited spectroscopy on the former. 

A more robust approach would be the derivation of two-point 
correlation statistics of the star forming galaxies, which could di- 
rectly test for clustering tendencies in the CDFS field compared to 
S 1 1 . Such analyses are beyond the scope of this paper, but will be 
addressed in a future work. 



6 SUMMARY AND CONCLUSIONS 

In this paper we report the results of a survey for Ha emitting galax- 
ies at z ~ 0.24. We used two fields from the Wide Field Imager 
Lyman Alpha Search (WFILAS). It consists of imaging in three 
narrowband filters (FWHM = 70 A), an encompassing intermediate 
band filter (FWHM = 220 A), supplemented with broadband B and 
R. The narrowband filters cover a redshift range of 0.23 < z < 0.26 
for Ha galaxies. These galaxies were selected by having an excess 
flux in one of the narrowband over to the other two, while also being 
detected in the intermediate and broadband R filters. This yielded a 
total of 707 candidate emission line galaxies (after the removal of 
stellar contaminants) for both fields. 

Of the 372 and 335 candidates, we observed 301 and 255 
through spectroscopic follow-up for the CDFS and Sll fields, re- 
spectively. We have identified emission in 189 and 117 candidates 
and confirmed that around half of these galaxies are Ha at z ~ 0.24. 
A significant number of galaxies were also found at z ~ 0.21 by 
means of their [S ii] emission. Other galaxies found were [On] and 
H/3/[Oiii] emitters at z ~ 1.2 and z ~ 0.6-0.7, respectively. Through 
use of the spectroscopy, we refined our colour selection to account 
for galaxies with a single emission line, leading to a measure of the 
fraction of Ha galaxies as a function of narrowband flux in both of 
these regions of the sky. We also used the spectroscopy to deter- 
mine a generic extinction correction using the Balmer decrement. 

We have determined the Ha luminosity function at z ~ 0.24 
separately for both of our fields after correcting for imaging and 
spectroscopic incompleteness, extinction and contamination from 
interlopers. We find small differences in their slope and turn-over 
luminosity while their normalisations were the same. When com- 
pared to recent Ha surveys, there is remarkable agreement between 
the luminosity function o f our CDFS field with th at one the Fabry- 
Perot imaging survey of lHippeleineTZI ( |2003|) . Differences be- 
tween our fields were of the order expected by cosmic variance but 
less than the scatter between the Ha luminosity functions of recent 
surveys. We surmise that while cosmic variance is a major contribu- 
tor to this scatter, it is differences in methodology between surveys 
(mainly differences in selection criteria) that dominate discrepan- 
cies between Ha luminosity functions and its related observables 
at z ~ 0.24. A survey that covers 10 - 20x the volume of one of 
our fields is re quired to get the unc ertainty due to cosmic variance 
to the levels of lCallego et"ai] ( Il995l) . 

We estimated the star formation density for both our fields 
to be logp = -1.93+°'^^ and -2.24+°[{ (p in Moyr"') for the 
CDFS and Sll fields, respectively, down to our survey limit of 
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logF,i„e = -16.0 (F,i„e in ergs ' cm -) or logLune = 40.6 {L,„,^ in 
ergs"'). Tliese values are comparable to other surveys at this red- 
shift when calculated to the same flux limit. Correcting for AGN 
would decrease these values by 0.02 to 0.04 depending on exactly 
how much of the Ha flux is contributed by the active nucleus rather 
than by normal star formation. 

Furthermore, we determined the star formation density in the 
hypothetical case where [S n] emitters atz ~ 0.21 were classified as 
Ha to illustrate the problems associated with solely relying colour 
selections. The star formation density logp of the SU field does 
not change by much (+0.02). On the other hand, the star formation 
density in the CDFS increases by 0.12, due to the large number of 
foreground [S ii] galaxies at z ~ 0.21. 

We explored the amount of star formation with respect to the 
local environment and found that the star formation rates were typ- 
ical for the field galaxy densities probed, in agreement with the re- 
sults of previous work. However, we also found tentative evidence 
of an increase in star formation rate per galaxy with increasing den- 
sity of the star forming galaxies. This supports scenarios where 
merger events are triggers for enhanced star formation, provided it 
can be demonstrated to be occurring on the smallest scales. We ex- 
plored this trend by examining the spatial distribution of our fields 
individually and found that it was largely attributable to one field. 
A formal study of the clustering statistics of this field is required to 
confirm this and will be the subject of a future study. 
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